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Abstract. In this paper, we generalize the parametric A-VaR method from 
portfolios with normally distributed risk factors to portfolios with elliptically 
distributed ones. We treat both the expected shortfall and the Value-at-Risk 
of such portfolios. Special attention is given to the particular case of a multi- 
variate t-distribution. 
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1. Introduction 

The original RiskMetrics methodology for estimating VaR was based on para- 
metric methods, and used the multi-variate normal distribution. This approach 
works well so-called linear portfolios, that is, those portfolios whose aggregate re- 
turn is, to a good approximation, a linear function of the returns of the individual 
assets which make up the portfolio, and in situations where the latter can be as- 
sumed to be jointly normally distributed. For other portfolios, like portfolios of 
derivatives depending non-linearly on the return of the underlying, or portfolios 
of non- normally distributed assets, one generally turns to Monte Carlo methods 
to estimate the VaR. Monte Carlo methodology has the obvious advantage of be- 
ing almost universally applicable, but has the disadvantage of being much slower 
than comparable parametric methods, when the latter are available. This is an 
issue in situations demanding for real-time evaluation of financial risk. For non- 
linear portfolios, practitioners, as an alternative to Monte Carlo, use A-normal 
VaR methodology, in which the portfolio return is linearly approximated, and an 
assumption of normality is made. Such methods present us with a trade-off between 
accuracy and speed, in the sense that they are much faster than Monte-carlo, but 
much less accurate unless the linear approximation is quite good, and the normal- 
ity hypothesis holds well. In case the linear approximation is of poor quality, or 
inherently instable (as is for example the case when the portfolio A is close to 0), 
one turns to a higher order approximations, for example the quadratic one, while 
keeping the normality assumption. This leads to the so-called F — A VaR, which 
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can be evaluated using Mont Carlo, but for which also a number of semi-parametric 
methods have been developed (see for example Duffy and Pan 2 , and others. 

An obvious first generalization is to keep the linearity assumption, but replace 
the normal distribution by some other family of multi-variate distributions. This 
is the subject of the present paper. As an alternative to the normal hypothesis we 
will assume that the log returns of the portfolio's constituents have a multivariate 
elliptic distribution in one of the elliptic classes N{^,T,, (p), cf. section 1 below 
for the precise definition. These have the advantage that, like normal distributions, 
their dependence structure is completely determined by their mean and variance, 
once a choice is made for 0. Particular examples are the normal distributions and 
the multi- variant Student-i distributions. The latter are an obvious first choice 
since they posses heavy tails. 

Glasserman, Heidelberger and Shahabuddin present a method to compute 
even F — A VaR using a semi-parametric method based on the Fourier transform, 
but their methodology seems to be restricted to t-distributions. See also Lopez and 
Walter [S] and references therein for further applications of the t distribution to 
VaR. Note that one shortcoming of the multivariate t-distribution is that all the 
marginal distributions must have the same degrees of freedom, which implies that 
every risk factor has equally heavy tails. 

In a sequel paper we intend to extend the present analysis to include quadrati- 
cally non-linear portfolios. We note that in a companion paper we give very precise 
analytic VaR estimates for the VaR of such quadratical portfolios, in case the un- 
derlying risk-factors follow another class of non-Gaussian distributions, namely the 
multivariate generalized normal distributions. 

The paper is organized, as follows: In section 2, we will analyze the VaR of a 
linear portfolio with elliptically distributed risk-factors, paying special attention to 
the case of a multi-variate Student distribution. Used in conjunction with a first 
order Taylor approximation of a portfolio's Profit &: Loss function, this will give rise 
to the notions of Delta- Elliptic VaR and Delta-Student VaR, in analogy with the 
familiar Delta-Normal VaR. We show, for example, how to reduce the computation 
of the Delta-Student VaR to finding the zero's of a certain special function. In 
section 3 we show how to extend our procedure to mixtures of elliptic distributions. 
Section 4 treats the expected shortfall for general elliptic linear portfolios and for 
the special case of Student ones. Finally, in section 5 we discuss some potential 
application areas. 

2. Linear Portfolio VaR with elliptic distributions 

In this section we perform a parametric analysis that relies on the assumption 
that the pricing function of the portfolio is linear in the risk factors. Note that 
parametric methods provide very fast answers which are, however, only as accurate 
as the underlying linearity assumption. 

We will use the following notational conventions for the action of matrices on 
vectors: single letters x,y, - ■ ■ will denote row vectors {xi, • • • , a;„), (yi, • • • j/n). 
The corresponding column vectors will be denoted by a:*,?/*, the * standing more 
generally for taking the transpose of any matrix. Matrices A = {Aij)ij, B , etc. 
will be multiplied in the usual way. In particular, A will act on vectors by left- 
multiplication on column vectors, Ay*' , and by right multiplication on row vectors, 
xA] X ■ X = XX* = Xi -\- ■ ■ ■ -\- x'^ will stand for the Euclidean inner product. 

A portfolio with time-t value 11 (i) is called linear if its profit and loss An(i) — 
n(<) —11(0) over a time window, [0 t] is a linear function of the returns Xi{t), . . . , V„ (t) 
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of its constituents over the same time period: 

An(t) = SiXi + 82X2 + ... + SnXn 

This will for instance be the case for ordinary portfolios of common stock, if we 
use percentage returns, and will also hold to good approximation with log-returns, 
provided the time window [0,t] is small. We will drop the time t from our notations, 
since it will be kept fixed, and simply write Xj,AIl, etc. We also put 

X = {Xi, ■ ■ ■ iXn), 

so that An = (5 - X = dX\ 

We now assume that the Xj are elliptically distributed with mean jj, and corre- 
lation matrix S = AA*: 

(Xi,...,X„)~7V(m,E,0). 
This means that the pdf of X is of the form 

fx{x) = \n-^g{{x-^i)Y.-\x-^y), 

where |S| stands for the determinant of S, and where g : K>o — > is such that the 
Fourier transform of g{\x\'^), as a generalized function on R", is equal to 
Assuming that g is continuous, and non-zero everywhere, the Value at Risk at a 
confidence level of 1 — a is given by solution of the following equation: 

(1) Prob {AU(t) < -VaRa} = a 

Here we follow the usual convention of recording portfolio losses by negative num- 
bers, but stating the Valuc-at-Risk as a positive quantity of money. 

In terms of our elliptic distribution parameters we have to solve the following 
equation: 



a = |Ef " / g{{x - fi)^-'{x - iiy)dx. 

Changing variables to y = {x — n)A~^ , dy =\ A \ dx , where E = A* A is a 
Cholesky decomposition of A, this becomes 

9{\y\^)dy. 

{6A-y<-6-iJ.-VaRc} 

Let i? be a rotation which sends SA to 0, 0). Changing variables once more 

to y = zR, we obtain the equation 

a= [ gi\z\')dz. 

J{\5A\zi<-5-n-VaRc} 

If we write that = -|- |.2'p with z' G M""^ then we have shown that : 

a = Pioh{6-X <-VaRa}= / [/ g{zl + \z'\^)dzi]dz' . 

Next, by using spherical variables z' = with ^ e Sn-2 , dz' = r"'~^da{^)dr, we 
see that we have to solve for VaRa in the equation 

-5n* -VaRg 



p+oo p 

a = \Sn-2\j^ r^-\j 



(2) a = \Sn-2\ / r"-^ / g{zi + r')dz^ 



-[TAT 



dr, 



"^One uses as a parameter for the class of elliptic distributions, since it is always well-defined 
as a continuous function: 0(|$P) is simply the characteristic function of an X ~ N[Q, Id, ij>)- Note, 
however, that in applications we'd rather know g 
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|iS'„„2| being the surface measure of the unit-sphere in 



l'S'n-2 



r(^) 



We now introduce the function 



4-C30 

dzi 







2 

r(iHI) 



f+OQ 

(3) = TTTTjirTT / {u ~ zl)"^ g{u)dudzi 



where for the second Hne we changed variables u = j-^ + z^. and replaced zi by 
— zi. We then have proved the following result: 

Theorem 2.1. Suppose that the portfolio's Profit & Loss function over the time 
window of interest is, to good approximation, given by AH — 5iXi + 82X2 + ••■ + 
dnXn, with constant portfolio weights 6j. Suppose moreover that the random vector 
X = (^1, • • • ,Xn) of underlying risk factors follows a continuous elliptic distribu- 
tion, with probability density given by fx{x) — |S| g[(x — ^)T,~^ (x — ^Y) where /x 
is the vector mean and E is the variance-covariance matrix, and where we suppose 
that (?(s^) is integrable over R, continuous and nowhere 0. Then the portfolio's 
Delta-elliptic VaR VaRa at confidence 1 — a is given by 



VaR^ = -S-^i + g9 „ • VST^, 
where Qa — n ^he, unique, solution of the transcendental equation 

a = G{qa,n)- 

Remark 2.2. Note that \dA\ has a clear financial interpretation, since 



(4) |M| = V5-S-(5*, 

which is simply the portfolio's volatility, or the square of its variance. 

Remark 2.3. In short-term Risk Management one can usually assume that /z ~ 0. 
In that case, we have that 



VaR^ = VS^ ■ 

which is completely analogous to the result for linear portfolios with normally dis- 
tributed risk factors, except that, for example for a — 0.05, the normal quantile 
at 5%, which is approximately 1.65, is now replaced by the g-dependent constant 
la 05- "^^^ latter will have to be computed numerically, for the different g's one 
would like to use. 

Remark 2.4. One can in fact do the integral over zi in 0: by Fubini, 

POO 

(5) G{s) = / K{s,u)g{u)du, 



where the kernel K is given by: 

K{s,u) = -|5'„_2| / (u-zi) 2 dzi 



1 3 1 

-|S'„_2| y (u~y)'^y-^dy 

-:\Sn-2\ x^{u-xyidx. 
4 Jo 
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At this stage we can the foUowing integral from Gradshteyn and Ryzhik |31 ■ 



-dx = — 2Fi{iy, fi; n + 



provided Re /i > and |arg + < tt; cf. 4 , formula 3.194(1). It follows that 

. . 7r~2~ Ti-i fln—ln+l 

(6) K{s,u) p(-»+i) ("~ ^) ' 2i^i (^2' ~2~' ^"'""^"^ 

However, we shall see in the example of the multi-variate t-distribution which we 
will treat next, that it can be easier to work directly with the double integral version 
© instead of with (0), ©. 

2.1. The case of t-student Distributions. We now consider in detail the case 
where our elliptic distribution is a multivariate Student-i. We will, unsurprisingly, 
call the corresponding VaR the Delta-Student VaR, generalizing the familiar ter- 
minology of Delta Normal VaR. 

In the case of multi-variate t-student distributions, the portfolio probability den- 
sity function is given by: 

X G M" and v > 2. Hence g is given by 

g{s)^C{v,n){l + s/v)-'-^ , s > 0, 

Where we have put 

C{v, n) 



r(^) 



Using this g in ©, we find that 



(7) G{s)^'^\Sn-2\C{v,n) I{zi)dz,, 

J S 

where we have put 

(u - zf)^{iy + uy~^du. 

-1 

The function I{zi) can be evaluated with the help of another one of the integrals 
inm: 

Lemma 2.5. (Cf. H] ,page 314-) If\arg{^)\ < n, and Re{vi) > Re{^) > , then 

/•+00 

(9) / ix-w)^'~\p + x)-'''dx = iw + pf-'''B{iyi-^i,fi), 

J W 

with B{a, (3) the Euler Beta function: 

r(Q; + P) 

Using formula Q with i^i — , /i — ^3^,/3 = ly, and w ~ zf , and therefore, 
II — i^i = —^-^ and — + vi — ii^, we find that 

(10) /(zi) = (z? + .)''^S(^,^) 
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We have not finished yet, since we still have to integrate over Zi in JT)). We 
therefore have to evaluate 



(11) Jis,iy) = I {zl + v) '^^'^dzi 
Changing variable in this integral according to u = z^, we find that 

1 r°° 1 -^i±ii 

(12) J{s,iy)^- / u^2(^u + iy) 2 du 

For the latter integral, we will use another formula fi'om 

Lemma 2.6. (cf. 4^, formula 3.194(2)). If\arg{^)\ <tt, and Re{vi) > Re{p) >0 
, then 

(13) / x'^-^l + Pxy^dx^ ^ ^F^{i^uiyi-fi;iyi-^i+l;-^). 

Ju Vi- fj. fj-U 

Here 2-Fi(a; 7; tu) is the hypergeometric function. 

In our case, 1^1 — ii^, fi — ^, i/i — jj, = ^,[3 — v^^ , and u = s^. If we replace in 
H12|l . we will obtain the following expression. 

(») V^n(i±i^,fa + ^.-i. 

Recalling |(7J), we find, after a small computation, that in the Student-t case, 
G{s)^G'^{s) = i^^|5„_2|C(^.,n).--2Fi(i±^,^;l + ^;-il) 



vJ^ \s'^) V{^) ' n 2 ' 2' 2' 



V 



(15) 

> 2 ^ 

Hence we have proved the following result on Delta-Student VaR: 

Theorem 2.7. Assuming that AH ~ 5iXi + 82X2 + ... + 5nXn with a multivari- 
ate Student-t random vector (Xi, X2, Xn) with vector mean 11 , and variance- 
covariance matrix S, the linear Value-at-Risk at confidence 1 — a is given by the 
following formula 

where now s — ql^ is the unique positive solution of the transcendental equation 
with Gi defined by U5\) . 

Remark 2.8. Note that gj, ^ does not depend of n. 

Hypergeometric 2^1 's have been extensively studies, and numerical software for 
their evaluation is available in Maple and in Mathematica. 



2.2. Some Numerical Result of Delta student VaR coefficient go!,y In the 
following table, we estimate only the positive solution of G{s) — a for some u 
Given, with the help of Mathematica 4 Software. 



V 


2 


3 


4 


5 


6 


7 


8 


9 


'Zo.oi,!/ 


6.96456 


4.54056 


3.74695 


3.36493 


3.14267 


2.99795 


2.89646 


2.8214 


90.025,1/ 


4.3026 


3.18244 


2.77644 


2.57058 


2.44691 


2.36462 


2.3060 


2.26216 


90.05,1/ 


2.91999 


2.35336 


2.13185 


2.01505 


1.94318 


1.89458 


1.85955 


1.81246 
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V 


10 


100 


200 


250 


275 


300 


400 


1000 


90. 01,1/ 


2.76377 


2.36422 


2.34135 


2.34514 


2.33998 


2.33884 


2.33571 


2.33008 


90.025,1/ 


2.22814 


1.98397 


1.97189 


1.96949 


1.96862 


1.9679 


1.96591 


1.96234 


90.05,1/ 


1.66023 


1.66023 


1.65251 


1.65097 


1.65041 


1.64995 


1.64867 


1.64638 



Remark 2.9. Note that, we obtain practically the same result as in the case of 
normal distribution, when the degree of freedom of our t-student is sufficiently high 
(y near 300), as it of course should. 



3. Linear VaR with mixtures of elliptic Distributions 



Mixture distributions can be used to model situations where the data can be 
viewed as arising from two or more distinct classes of populations; see also 0. For 
example, in the context of Risk Management, if we divide trading days into two sets, 
quiet days and hectic days, a mixture model will be based on the fact that returns 
are moderate on quiet days, but can be unusually large or small on hectic days. 
Practical applications of mixture models to compute VaR can be found in Zangari 
(1996), who uses a mixture normal to incorporate fat tails in VaR estimation. 
Here we sketch how to generalize the preceding section to the situation where the 
joint log-returns follow a mixture of elliptic distributions, that is, a convex linear 
combination of elliptic distributions. 

Definition 3.1. We say that {X\, ...,X„) has a joint distribution that is the mix- 
ture of m elliptic distributions N{iij, Sj, (pj)'^, with weights {Pj} (j=l,..,m ; /3j > 
; X^jLi = 1); if it^ cumulative distribution function can be written as 

m 

Fxi,...,X„ixi, ...,Xn) = ^f3jFj{xi, ...,Xn) 

with Fj{xi, Xn) the cdf of N{iij, Sj, (f>j). 

Remark 3.2. In practice, one would usually limit oneself to m = 2, due to esti- 
mation and identification problems; see 6_. 

We will suppose that all our elliptic distributions N{pj, Sj, (f>j) admit a pdf : 



(16) 



f,ix) = \E,\-'/^g,{{x - - M,)*)- 



The pdf of the mixture will then simp;y be X^^i f^jfji^ 



Let 



Ej — A* Aj 



be a Cholesky decomposition of Ej . Since integration is a linear operation, we now 
have to solve 

-6-iJ,j -VaRa 

(17) a=\S„^2\Y.l3,\^,\-'/' "^^ 



r+oo 












Jo 





dr 



to obtain VaRa- This leads to the following theorem: 



'or N{fj,j, Hj,gj) if we parameterize elliptical distributions using g instead of (/> 
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Theorem 3.3. Let AH = SiXi + . . . + (5„X„ with {Xi, . . . , Xn) is a mixture of 
elliptic distributions, with density 

m 

fix) = ^/3,|I],ri/2g,((x - - fi.Y) 

where is the vector mean, and Y.j the variance- covariance matrix of the j-th 
component of the mixture. We suppose that each gj is integrable function over R, 
and that the gj never vanish jointly in a point of M™. Then the value-at-Risk, 
or Delta mixture-elliptic VaR, at confidence 1 — a is given as the solution of the 
transcendental equation 

where Gj is defined by ^ with g — gj. Here S — (Si, ... , (5„). 

Remark 3.4. In the case of a mixture of m elliptic distributions the VaR will 
not depend any more in a simple way on the total portfolio mean and variance- 
covariance. This is unfortunate, but already the case for a mixture of normal 
distributions. 

Remark 3.5. One might, in certain situations, try to model with a mixture of 
elliptic distributions which all have the same variance-covariance and the same 
mean, and obtain for example a mixture of different tail behaviors by playing with 
the gj's. In that case the VaR again simplifies to: VaRa = —S ■ fJ, + Qa ■ V 
with Qa now the unique positive solution to 

m 

a = ^l3jGj{qa). 

The preceding can immediately be specialized to a mixture of Student t-distributions: 
the details will be left to the reader. 



4. Expected Shortfall for elliptic distributions 



Expected shortfall is a sub-additive risk statistic that describes how large losses 
are on average when they exceed the VaR level. Expected shortfall will therefore 
give an indication of the size of extreme losses when the VaR threshold is breached. 
We will evaluate the expected shortfall for a linear portfolio under the hypothe- 
sis of elliptically distributed risk factors. Mathematically, the expected shortfall 
associated with a given VaR is defined as: 

Expected Shortfall = E(-An| - AH > VaR), 

see for example (101. Assuming again a multivariate elliptic probability density 
f{x) = |S]|~ g{{x ~ ii)Y.~^ [x — nY) , the Expected Shortfall at confidence level 1 — a 
is given by 

-ESa = E(An I An < -VaRo,) 
= iE(An- l{An<-yafl„}) 



1 

a 

IS 



(5a;* f{x) dx 

{Sx*<-VaRc,} 
-1/2 

5x^ g{{x — ii)'S^^{x — 

{Sxt<-VaR^} 



a 
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Let E = A* yl, as before. Doing the same linear changes of variables as in section 2, 
we arrive at: 

1 

a 
1 



^ESa = -/ {\8A\zx^8 ■ [i) g{\\zf)dz 



|(5A|zi (jr(||z||^) dz + (5 • /i. 

The final integral on the right hand side can be treated as before, by writing ||2;||^ = 
zj + ||z II and introducing spherical coordinates z = r^, ^ G Sn-2, leading to: 



oo 



\5A\ zi g{zi + r^)dzi dr 



We now first change zi into — zi, and then introduce u — z^ + , as before. If we 
recall that, by theorem l2.1l 



g _ S ■ ^ + VaRa 

9a, 



|M| 

then, simply writing for „, we arrive at: 

ESc, = -S-fi+\SA\^-^^- f f zi{u - zf)"^ g{u) du dzi 



oo 



\Q I r°° 1 71-1 
= _5.^+|M|^^./ -{u-ql)^ giu)du, 

since 

^ Z,{u-zl)^--' = -^^{u-qlY-^ . 

After substituting the formula for |S'n-2| and using the functional equation for the 
F- function, r(x + 1) = xT{x), we arrive at the following result: 

Theorem 4.1. Suppose that the portfolio is linear in the risk-factors X = {Xi, • • • , Xn): 
An = S-X and that X ^ N{iJ,,T,,(j)), withpdf f{x) = g{{x - {x - ^if) . 

If we write qa =, then the expected Shortfall at level a is given by : 



(19) ESo,^-6-t,+ \6J:6r'- {u-iqlnY) ' 9{u)du 



4.1. Application: Student Expected Shortfall. In the case of multi-variate 

t-student distributions we have that g{u) = C{v,n){\ -\-u/v) ^ , with C{v,n) 
given in section 2. Let us momentarily write q for j^. We can evaluate the integral 
in 1)19(1 using lemma VU\ as follows: 



du 



If we pose that : 

_ ^ ^ \ 



r(f) ^^^"'"-^ ^ > 



a ■ 

After substitution in (|19|l . we find, after some computations, the following result: 
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Theorem 4.2. The Expected Shortfall at confidence level 1 — a for a multi-variate 
Student-distributed linear portfolio S ■ X, with 

r(^) , (x-m)*s-i(x-m)\"^'^^ 



X 



1 + 



^5 given by: 



a-V^ r(|) 
-5- iJ, + eSa,u-\Si:S*\^/'^ 



|JE(5|V2 



The Expected Shortfall for a linear Student portfolio is therefore given by a com- 
pletely explicit formula, once the VaR is known. Observe that, as for the VaR, the 
only dependence on the portfolio dimension is through the portfolio mean 6 ■ /U and 
the portfolio variance 6T,6*. 

With the lu^lp of Alatlab wc olitaiii the tliv Mowing tabic of A'ahu^s of f ,,.,/: 



z/ 


2 


3 


4 


5 


6 


7 


8 


9 


eso.01,1/ 


5.5722 


5.9309 


5.7879 


5.4555 


5.0799 


4.7160 


4.3819 


4.0818 


eSo.025,!^ 


8.6113 


7.6777 


6.8216 


6.0676 


5.4326 


4.9032 


4.4601 


4.0862 


eso.05,1/ 


11.7123 


9.0750 


7.4966 


6.3797 


5.5457 


4.9007 


4.3880 


3.9711 



V 


10 


100 


200 


250 


eso.01,1/ 


3.8135 


0.5157 


0.2644 


0.2086 


eso. 025,1/ 


3.7675 


0.4577 


0.2313 


0.1854 




3.6257 


0.4073 


0.2050 


0.1642 



4.2. ES for Elliptic distribution Mixtures. The preceding results can, as be- 
fore, easily be generalized to mixtures of elliptic distributions. Details will be left 
to the reader. 

5. Some Areas of Applications 

In this section, we survey some areas for applications of linear portfolios that 
exist in the financial literature. We will discuss 4 examples: 

• Delta-approximation of a derivatives portfolio. 

• Linear approximation of an equity portfolio 

• Businesses as portfolios of business units 

• Incremental VaR 

5.1. Delta Approximation of a Portfolio. Suppose that wc arc holding a port- 
folio of derivatives depending on n underlying assets X^^^, X^^^ X("With el- 
liptically distributed log-returns r^ (over some fixed time- window) . The portfolio's 
present value Y will in general be some complicated non- linear function of the X^'s. 
To obtain a first approximation of its VaR, we simply approximate the present Value 
V of the position using a first order Taylor expansion: 

" dV 

V{X + AX) « V(X) + ^ ^AX,. 
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^From this, we can then approximate the profit/loss function as 

n 

AV = V{X + AX) - V{X) » ^ 5,r« 

i=l 

where we put r — r^"^) and S — (^i, (5„) with Si = X''^ g^^. The entries 

of the S vector are called the "delta equivalents " for the position, and they can 
be interpreted as the sensitivities of the position with respect to changes in each of 
the risk factors. For more details see where a multi-variate normal distribution 
for the Ti's is assumed. The discussion there generalizes straightforwardly to the 
elliptic case, where the present paper's results can be used. 

5.2. Portfolios of Equities. A special case of the preceding is that of an equity 
portfolio, build of stock Si, ■ ■ ■ , Sn with log-returns ri(t), •, r„(t). In this case, 

n 

n(0-n(0) = ^u;,^,(0)(5,(t)/5,(0)-l) 

n 
i=l 

where this approximation will be good if the (t) are small. 

5.3. Businesses as Linear Portfolios of Business Units. An interesting way 
of looking upon an big enterprize, e.g. a multi-national or a big financial institution, 
is by considering it as a sum of its individual business units, cf. Dowd |7j. If Xj, is 
the the variation of price or of profitability of business unit j in one period , then 
the variation of price of the agglomerate in the same period will be 

AU = Xi + ---+Xn. 

The entire institution is therefore modelled by a linear portfolio, with i5 = (1, 1, . . . , 1) 
to which the results of this paper can be applied, if we model the vector of indi- 
vidual price variations by a multi-variate elliptic distribution. VaR, incremental 
VaR (see below) and Expected Shortfall will be relevant here. For more details see 
Dowd [3, chapter XI . 

5.4. Incremental VaR. Incremental VaR is defined in [H] as the statistic that 
provides information regarding the sensitivity of VaR to changes in the portfolio 
holdings. It therefore gives an estimation of the change in VaR resulting from a 
risk management decision. Results from 10 for incremental VaR with normally 
distributed risk-factors generalize straightforwardly to elliptically distributed ones: 
if we denote by IVaRi the incremental VaR for each position in the portfolio, with 
9i the percentage change in size of each position, then the change in VaR will be 
given by 

AVaR = e jVaR, 
By using the definition of IVaRi as in 6 (2001), we have that 

(20) iVaR,=u;,^^ 

with uJi is the amount of money invested in instrument i. In the case of an equity 
portfolio in the elliptically distributed assets, we have seen that, assuming /i = 0, 



VaR^ - -q9^„V^, 
We can then calculate IVaRi for the i-th constituent of portfolio as 

I VaRi = = LOi^i 

duJi 
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with 

= -9 

The vector 7 can be interpreted as a gradient of sensitivities of VaR with respect 
to the risk factors. This is the same as in except of course that the quantile has 
changed from the normal one to the one associated to g. 

6. CONCLUSION 



In this paper we have shown how to reduce the estimation of Value-at-Risk for 
hnear elhptic portfoUos to the evaluation of one dimensional integrals which, for 
the special case of a Student distribution, can be explicitly evaluated in terms 
of a hypergeometric function. We indicated how to extend these to the case of 
mixtures of elliptic distributions. We have also given a similar, but simpler, integral 
formula for the expected shortfall of such portfolios which, again, can be completely 
evaluated in the Student case. We finally surveyed some potential application areas. 
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